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Abstract. The Cepheid mass discrepancy, the difference between mass estimates 
from stellar pulsation and stellar evolution models, is a long standing challenge for the 
understanding of stellar astrophysics. We discuss the current state of the mass discrep- 
ancy and test the role of pulsation-driven mass loss in state-of-the-art stellar evolution 
calculations of Large Magellanic Cepheids in resolving it. We find that Cepheid mass 
loss is a significant contributor to the mass discrepancy, but it is not clear if the metal- 
licity dependence of Cepheid mass loss is consistent with the measured metallicity de- 
pendence of the mass discrepancy. 



1. Introduction 

Classical Cepheids are ideal laboratories for understanding stellar astrophysics and evo- 
lution. Observations of Cepheid pulsation probes the interior structures of these stars 
constraining stellar evolution theory, while their luminosities logL/L© > 3 allows us 
to observe Cepheids in distant galaxies with differing metallicities. Furthermore, the 
tight correlation between the pulsation period and stellar luminosity also makes them 
powerful standard candles. However, for all their importance for stellar astrophysics 
and cosmology, the masses of Cepheids are still not understood. 

Cepheid masses are determined in three ways: from stellar evolution models, stel- 
lar pulsation models, and from measurements of binary systems with a Cepheid com- 
ponent. For example, iPietrzynski et all (l2010h determined the mass of the only known 
fundamental-mode Cepheid to be in an eclipsing binary. Other dynamical masses have 
been determin ed for Cepheids in othe r types of binary systems, but not to the same 
precision (e.g. lEvans et alllloOSL l2009h . Masses from stellar evolution models are de- 
termined by fitting evolutionary tracks to the observed effective temperature and lu- 
minosity, while masses from stellar pulsation models are determined by fitting pulsa- 
tion properties. The resulting paradox is that these latter two methods pre dict d ifferent 
masses, this is historical problem is called the Cepheid mass discrepancy (ICoxl|l9 80). 

Currently, the Cepheid mass discrepan cy has been m easured for Galactic, and 
Large and Small Magellanic Cloud Cepheids. ' Keller & Wo od (2006) and Keller (2008]) 
determined that the mass discrepancy ranges from about 15-25%, where the mass dis- 
crepancy is defined as the relative difference between evolutionary and pulsation predic- 
tions. Furthermore, they find that the mass discrepancy increases decreasing metallicity. 
It should be noted that stellar evolution masses in these works are computed assuming 
no extra mixing processes such as rotational mixing or convective core overshooting. 
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Bono et all ( 20061) reviewed possible solutions to the Cepheid mass discrepancy, 



where the two most likely candidates ai^e mass loss during the Cepheid stage of evolu- 
tion and convective core overshooting in the main sequence progenitors of Cepheids. 
Convective core overshooting is defined as the distance above the stellar core that con- 
vective eddies penetrate and mix over an evolutionary timescale. This distance is typ- 
ically defined as A = ctcHp, where Hp is the pressure scale height and is a free 
parameter. Overshooting mixes additional hydrogen into the core of a main sequence 
progenitor, producing a more massive helium core during the Cepheid stage of evolu- 
tion. Therefore, by including overshooting one can fit observed luminosity of a Cepheid 
with a stellar evolution track assuming smaller masses than an evolution model with 
no overshooting, and thus brings evolution calculations into agreement with pulsation 
models. However, Keller (2008.) found that a mass discrepancy of 20% requires a value 
of Oc = 0.8. This v alue of Or is inconsistent w ith observations of other stars such 



as ec lipsing binaries (San dberg Lacv et al 



2010l) . and massive B-type stars ( Brott et al. 



2010) , j6 Cephei stars (iLovekin & Goupil 

2011) . 



Mass loss is another possible solution; 



Neilson & Lested(l2008L 120091) derived an 



analytic theory describing pulsation-driven mass loss in Classical Cepheids, and found 
that mass-loss rates could be enhanced by more than three orders-of-magnitude rel- 
ative to a radiati on driven wind. Observat ional evidence for pulsation-driven mass 
loss is growing; iNeilson et all (I2009L I20T0I) measured mass-loss rates from infrared 
excess for LMC Ce pheids in the OGLE-II/II I surveys combined with the SAGE s urvey 
(iNgeow & Kanbuii f2008; Ngeo w et al.ll2009l) . JVIarengo et al. (2010) and Barmby etal] 
(I2OI ll) determined mass-loss rate s, M = 10~^ to 10~^ M p, yr~', for Galactic Cepheids 
from infrared observations, while iMatthews et all (1201 ih measured a mass-loss rate of 
10~^ Mq yr~' for 6 Cephei, based on 21 cm observations. This evidence is pointing 
towards enhanced mass loss in Cepheids and is thus worth testing if Cepheid mass loss 

resolves the mass discrepancy. 

In this work, we compute stellar evolution models including the iNeilson & Lester 

(120081) Cepheid mass loss prescription wit h standard LMC m etallicity and assuming 
moderate convective core overshooting from Brott et al.l ( 201 ih . Oc = 0.335, to compare 
with the results from .Keller & Wood (2006). The models are computed through blue 
loop evolution and we determine the change of mass due to Cepheid mass loss. We 
test if this amount of overshooting and mass loss is consistent with the measured mass 
discrepancy. Furthermo re, we compare the r esults from LMC Cepheids with the results 
for Galactic metallicity dNeilson et al.ll2011 ). 



2. Method 



We compute stella r evolution mod els using the lYoon & Langerl (120051) version of the 
state-of-the-art He ger et all (l2000l) s tellar evolution code. I n this code, we include a 
prescription for Cepheid mass loss (^ Neilson & Lestetll2008 ') that computes pulsation- 
driven mass loss when the effective temperature and luminosity of the stellar model 
falls within the Cepheid instability strip as defined by iBono et ahl (l2000l) . The mass- 
loss theory predicts the mass-loss rate as a function of luminosity, mass, radius, pul- 
sation period and pulsation amplitudes. The pu lsation period is determined from the 



period-mass-radius relation ( Gieren et aL 19891). and the pu lsation amplitudes for lu- 



minosity and radius are from Klagy ivik & SzabadosI (l2009h . Details can be found in 
INeilson et all (I2OIII) . We note that these relations are derived for Galactic Cepheids, 
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Figure 1 . Stellar evolution tracks for models with Galactic metallicity, Z = 0.02 
for two scenarios. The first (left) is for models assuming no mass loss a nd Or = 0, 
and the second (right) is for stellar evolution models including the N eilson & Lested 
12008) Cepheid mass loss prescription and Uc - 0.335. The dashed lines outUne the 
boundaries of the Cepheid instability strip dBono et al.ll2000h . 



hence Galactic metallicity, and thus ignores any implicit metallicity dependence of the 
pulsation period and amplitudes. 



3. Results 



In lNeilson et al.l (120111) . we compute stellar evolution models with Galactic metallicity, 
Z = 0.02, and masses M = 4 to 9 Mq in steps of 1 Mq for four scenarios. The first 
scenario is no mass loss and no convective core overshooting, while remaining scenar- 
ios include pulsation-driven mass loss and convective core overshooting with values of 
Oc = 0,0.1, and 0.335. We show the stellar evolution tracks for models with no mass 
loss and zero convective core overshooting include along with models with pulsation- 
driven mass loss included and = 0.335 in Fig. [T] Stellar evolution tracks without 
mass loss and convective core overshooting tend to have wider blue loops than model 
with Cepheid mass loss and = 0.335 for the same initial mass, where the width is 
defined by the change of effective temperature across the blue loop. One might con- 
sider that the shape of the blue loops is an argument against C epheid mass loss, b ut, the 
question of why a blue loop occurs at all is still unanswered. IValle et all ( 20091) found 
that the blue loop structure is alfected by the assumed stellar m ass, composition, and 
amount of convective core overshooting. We refer the reader to iNeilson et all (1201 ll) 
for more discussion. 



We repeat the calculations from Neilson et all (1201 11) for s tellar evolution models 
with LMC metallicity, Z = 0.008. In the Yoon & Langed (l2005h stellar evolution code, 
the helium abundance is correlated with the metallicity, therefore LMC stellar evolution 
models also assume a smaller helium abundance than Galactic models. Evolution mod- 
els are computed assuming = 0.335 and initial masses, M = 4,5, 6, and 7 Mq, with 
evolution tracks shown in F i g. El a long with the estimated boundaries of the instability 
strip defined by iBono et all (l200Clh . These are the bound aries determined for Galactic 
Cepheids, however we demonstrated ( Neilson et al.ll20lil) that the total amount of mass 
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Figure 2. Stellar evolution tracks for models assuming the standard LMC metal- 
licity, Z = 0.008. The dashed lines denote the assumed boundaries of the Cepheid 
instability strip. 



lost during the blue loop evolution is largely insensitive to assumed location of the red 
edge of the instability strip. 

The mass-loss theory developed by lNeilson & Lesteil (120081 '). based on the lCastor et al, 

(| 19751 ) CAK radiative-driven wind theory, noted that Cepheid mass-loss rates explic- 
itly depend on metallicity, M oc (Z/0.02)^''^. Therefore, we would expect that LMC 
Cepheid mass-loss rates to be about 0.63x mass-loss rates for Galactic Cepheids, if 
all other properties are the same. Assuming the evolutionary timescale of the Cepheid 
blue loop is the same for both Galactic and LMC metallicities then we would expect 
the Cepheid mass discrepancy for LMC Cepheids to be about 3 - 4% smaller than that 
for Galactic Cepheids. Furthermore, since the LMC models have smaller helium abun- 
dances then the luminosity during blue loop evolution is also smaller. Therefore, the 
expected difference between Galac tic and LMC ma s s discr epancies would be > 3-4%, 
while the observed difference from lKeller & WoodI (l2006l) is about -3%. 

In Fig.[3l we show the predicted mass discrepancy for stellar evolution models with 
Galactic and LMC metallicities for the assumed stellar masse s in this work . We also 
plot the measured mass discrepancy for Galactic Cepheids from lKelleil ( 2008h . 17 ±5%, 
as a shaded box to account for the uncertainties. The mass-loss rates for LMC Cepheids 
tend to be smaller than for Galactic Cepheids leading to a smaller Cepheid mass dis- 
crepancy with decreasing metallicity. The difference between predicted Galactic and 
LMC mass discrepancy ranges from about zero at 4 Mq to 7% at 6 Mq. However, the 
difference is not a function of initial mass, at 7 Mq the difference is about 2%. 

The results presented in Fig. |3] are not consistent with the metallicity dependence 
of the mass discrepancy determined by lKeller & Wood! (l2006h . However, the predicted 
mass discrepancy as a function of metallicity is different than naively expected, suggest- 
ing that differences in stellar properties due to assumed metallicity are acting to predict 
larger mass-loss rates at smaller metallicities. Another issue is the use of the period- 
pulsation amplitude relations calibrated using Galactic Cepheids; iBono et all (l2000l) 
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Figure 3. Predicted Cepheid mass discrepancy for Galactic (magenta circles) and 
LMC (blue squares) as function of initial stellar mass. The shaded region represents 
the measured Cepheid mass discrepancy of 17 ± 5% for Galactic Cepheids (Kellerl 
l20Q8ft . 



found that pulsation mode ls of LMC Ceph eids have larger amplitudes than similar mass 
Galactic Cepheids. From iNeilson & Les ter (2008), an increase of the brightness am- 
plitude by a factor of 2 leads to an increase of mass loss by about 5%, which leads 
to a significantly increased predicted mass discrepancy at LMC metallicity. This does 
not even include the effect of a larger radius amplitiide on the predicted mass-loss rates. 
This behavior is consistent w ith the resul ts of lNeilson & Lested (j2009t) . where mass-loss 
rates are computed for Bono et all ( 20001) theoretical models and found larger mass-loss 
rates for LMC composition relative to Galactic metallicity. 



4. Summary 

In this work, we coni puted stellar evolution models with L MC metallicity us i ng the 
Yoon & Langed (l2005l) stellar evolution code that includes the lNeilson & Lesteil (l2008l) 



prescription for Cepheid mass loss. From these models, we tested how much of the 
measured Cepheid mass discrepancy (.Keller & Wood.2006 ) can be accounted for by as- 
suming moderate convective core overshooting and pulsation-driven mass loss. We also 
comp ared these prediction s with mass discrepancy predictions for Galactic Cepheids 
from iNeilson et al] ( 201 ih to see if we a gree with the metallicity dependence of the 



mass discrepancy found by lKeller & Woo d (2006). 

The results are not precise enough to demonstrate that pulsation-driven mass loss 
plus moderate convective core overshooting, = 0.335, is the solution to the metal- 
licity dependence of the Cepheid mass discrepancy. However, by assuming the pul- 
sation amplitudes are the same for Galactic and LMC metallicities, we have underes- 
timated both the pulsation amplitudes and mass-loss rates. Thus, we underestimated 
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the contribution of Cepheid mass loss to the LMC Cepheid mass discrepancy, hence, 
we cannot rule out the possibility that pulsation-driven mass loss is consistent with the 
results of Kelle r & Wood (2006). This is encouraging since one would naively expect 
significantly smaller mass-loss rates for smaller metallicities, which we do not find. 
Pulsation-driven mass loss plus convective core overshooting provides a plausiable so- 
lution to the Cepheid mass discrepancy, but more precise studies are required to verify 
this hypothesis. 

Acknowledgments. HRN is grateful for funding from the Alexander von Hum- 
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